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ABSTRACT 

We have developed a new method for calculating common envelope (CE) events based on explicit 
consideration of the donor star's structural response to adiabatic mass loss. In contrast to existing 
CE prescriptions, which specify a priori the donor's remnant mass, we determine this quantity self- 
consistently and find it depends on binary and CE parameters. This aspect of our model is particularly 
important to realistic modeling for upper main sequence star donors without strongly degenerate cores 
(and hence without a clear corc/cnvelopc boundary). We illustrate the central features of our method 
by considering CE events involving 10 Mq donors on or before their red giant branch. For such 
donors, the remnant core mass can be as much as 30% larger than the star's He-core mass. Applied 
across a population of such binaries, our methodology results in a significantly broader remnant mass 
and final orbital separation distribution and a 20% increase in CE survival rates as compared to 
previous prescriptions for the CE phase. 

Subject headings: binaries: close ~ stars: evolution 



1. INTRODUCTION 

The common envelope (CE) phase of binary evolution 
has become a standard paradigm for the formation of 
short period systems containing a compact object since 
it was introduced by Paczynski (1976) (see reviews by 
Iben & Livio 1993; Taam & Sandquist 2000). While we 
lack quantitative criteria for determining whether the en- 
velope is successfully ejected in a CE event and a short- 
period binary system formed, prescriptions have been 
proposed to approximate the relation between the prop- 
erties of the post-CE binary and its pre-CE progenitor 
system. The major source of uncertainty in these pre- 
scriptions is parameterized by a factor, acE, denoting 
the efficiency by which orbital energy is converted to the 
kinetic energy of the ejected envelope (see Iben & Tu- 
tukov 1985) or a factor, 7, describing the specific angu- 
lar momentum of the ejected matter (see Nelemans et al. 
2000; Nelemans & Tout 2005). 

In recent years, significant effort has been devoted to 
understanding the early evolution of CE events in three 
dimensions (e.g. Ricker & Taam 2008) in order to elu- 
cidate the relevant hydrodynamic processes. However, 
a detailed description of the final phase, where the en- 
velope is fully ejected and the remnant binary system 
forms, is difficult to obtain via hydrodynamic simula- 
tions since resolving the spatial scales of the donor's core 
required prohibitively short time steps. The main goal of 
this letter is to describe a physically motivated yet com- 
putationally tractable methodology for determining the 
end-state of this final CE phase. 

The outcome of a CE event depends not only on the 
rate at which mass is ejected, but also on the ability of 
the donor's remnant layers to contract in response to the 
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rapid mass loss. The initial response of stars to adia- 
batic mass loss was carried out in the pioneering work 
of Hjellming (1989), where the central focus was on de- 
termining when dynamical mass transfer instabilities oc- 
cur. However the donor's response down to the removal 
of nearly the entire envelope (relevant to the final phase 
of CE evolution) was not considered. Here we report 
on our efforts to extend the work of Hjellming (1989) 
by determining the structural response of donor stars to 
adiabatic mass loss sequences that terminate deep within 
the donors' core. In §2, we briefly describe these compu- 
tations and their results. In §3, we apply these results 
to develop a new procedure for determining the outcome 
of the CE phase. We compare the results of an analy- 
sis for a binary population of a 10 Mq donor with those 
based on previous prescriptions for the properties of the 
remnant core in §4 and summarize in the last section. 

2. METHODOLOGY AND DONOR RESPONSE 

The structure and response of a donor star is calcu- 
lated using a code developed and described in Deloye 
et al. (2007) in the limit of adiabatic mass loss (i.e., 
we assume that in the CE, mass loss proceeds at a suf- 
ficiently high rate that energy transport between mass 
layers can be neglected). The equations governing the 
donor's structure are given by the condition of mechani- 
cal equilibrium subject to the constraint that the entropy 
of a given mass layer is fixed. The set of equations used 
to calculate the donor's structure is equivalent to the fol- 
lowing: 
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where p, T, P, r, and s arc the density, temperature, 
pressure, radius, and specific entropy for each mesh shell. 



and m is the mass interior to it. The initial entropy 
Si{m) and composition profile, X(m)^ is determined from 
a donor's initial structure. 

We apply standard boundary conditions except that 
in lieu of the effective temperature-luminosity relation, 
we enforce the fixed-entropy condition at the exterior 
boundary. To enforce mass loss on the donor's structure, 
the outer mesh point m equals the varying total donor 
mass, M2. 

In reality, the above set of equations is cast in terms of 
modified P, T, 771, and r variables and supplemented by 
a mesh spacing function that automatically adjusts the 
non-Lagrangian mesh as the evolution proceeds (Eggle- 
ton 1971). This recasting was adopted so that our equa- 
tions had the same form as the stellar evolution code 
used to calculate the initial models, providing for a more 
robust ability to reconverge these initial models. Along 
a mass loss sequence, the mass at each grid point varies. 
To calculate Si{m) and the composition, X(m) at each 
grid point, we interpolate the initial profiles using Akima 
splines. 

While we focus here on donors with zero-age main se- 
quence (MS) masses of M2,o ~ 10 M©, we calculated the 
adiabatic mass-loss response of donors with M2.0 in the 
range between 1.5-20 M© with the mass loss starting at 
evolutionary phases from the terminal age MS through 
the red giant and asymptotic giant branch phases. The 
models from which we initiated mass-loss sequences we 
calculated using the STARS code (Eggleton 1971, 1972, 
1973; updated as described in Pols et al. 1995) for a 
metallicity of Z = 0.02. We reconverged each input initial 
model (to an accuracy of ~ 10~^ in a norm of summed, 
squared relative residuals) in our own code with M2 held 
fixed. We then calculated the mass loss sequences by 
decrementing M2 by a small amount AM2, based on both 
M2 and the current surface entropy gradient, and con- 
verging the model at this new M2 value. The results of 
the full study will be presented in a follow-up paper. 

Our M2fl — 10 M0 mass loss sequences were calcu- 
lated with a numerical resolution of 599 grid points. The 
variation of stellar radius, i?2, with M2 for several of our 
sequences is shown for donors at the end of the MS, near 
the end of the Hcrtzsprung gap (HG) , and within the red 
giant branch (RGB) in Fig. 1. The marked difference in 
the initial response between the RGB model's mild ex- 
pansion and the other phase's rapid contraction is due 
to the transition from a radiative envelope in the less 
evolved donors to the development of a nearly flat en- 
tropy profile in the RGB stars' convective envelopes. As 
the donor moves up the RGB, the H— >He burning shell 
thins and increases in luminosity. This, in tandem with 
the resulting deepening of the convective envelope, pro- 
duces a steepening entropy profile in the core/envelope 
transition region. The steepening of the R2 response near 
the core/envelope boundary for more evolved models re- 
sults from this entropy gradient. Correspondingly, the 
more gradual contraction experienced during the later 
evolution of the less evolved donors results from a shal- 
lower entropy gradient throughout the donor's interior. 
The intersection of curves at m w 2.3 M© corresponds 
to the mass at the top of the H-burning shell. 

We emphasize that in the adiabatic limit, the R2 re- 
sponse to mass loss depends only on the donor's initial 
structure through its s-profile. Thus, in a CE event, the 
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Figure 1. The radius evolution of 10 M© donors undergoing 
mass loss for various initial evolutionary states. The dot-dash and 
dash-line curves (with radius at the onset of the rapid mass loss 
episode, R2,o = 10.8 and 163.3 R©) illustrate the adiabatie re- 
sponse of donors either near the end of the main sequence or near 
the end of the HG stage respectively. The five solid line curves 
(with _R2,o = 177.0, 216.2, 251.4, 287.2, and 318.7 R©) illustrate 
the response of donors at various stages on the RGB. 

i?2 vs. M2 evolution is independent of other external fac- 
tors (such as binary parameters) under the assumption 
that the binary companion does not affect the donor's 
structure. This allows us to calculate in advance a li- 
brary of R2~M2 tracks that form the essential input for 
an improved CE event outcome determination method- 
ology. 

3. A SELF-CONSISTENT METHOD FOR DETERMINING 
CE OUTCOMES 

During a CE phase, mass from the donor's envelope is 
stripped and then ejected from the binary at the expense 
of the orbital energy. Current prescriptions for the final 
outcome of this process, in terms of post-CE orbital sepa- 
ration, a/, take as an a priori input the donor's post-CE 
remnant mass, M2J. This mass is typically taken to be 
the mass of the donor's evolved core as determined via 
one of several criteria applied to isolated stellar evolu- 
tion models. For lower MS stars on their giant branches, 
this approximation is reasonable given the steep entropy 
profiles across their thin-shell burning regions. However, 
stars on the upper main sequence can have significantly 
shallower entropy gradients in the core/envelope bound- 
ary region. This leads to ambiguity in defining the core 
mass in these cases and produces a significant variation 
in predicted CE event outcomes (see, e.g., Tauris & Dewi 
2001). 

More importantly, in such cases, the physical motiva- 
tion for correlating core mass with M2J no longer holds. 
That is, the adiabatic mass loss sequence for lower MS 
donors have R2-M2 tracks that are similar to the track 
starting from the most-evolved donor with largest i?2,o 
in Fig. [3 where the steep drop-off in R2 near the core es- 
sentially guarantees the donor contracts within its Roche 
lobe with M2 nearly equal to the core mass. For cases 
with more gradual R2 evolution this is not guaranteed 
to be the case and the basic premise of using an a priori 
determined M2J becomes suspect. 

To address this issue, we have developed a new method 
for determining CE event outcomes based not on spec- 
ifying M2./, but on the R2 evolution. Specifically we 
assume that the CE event ends at the point the donor 



contracts within a critical radius, which we take as its 
Roche radius. To determine this point, both R2 and the 
donor's Roche lobe radius, Rl, are calculated as a func- 
tion of M2 during the CE event. Our mass loss R2-M2 
tracks provide the first necessary input. To determine 
the RL-M2 track during the CE event we assume the 
standard CE energy conservation (i.e., acE) prescription 
for determining aj applies along the entire course of the 
CE event. That is, the orbital separation, a during the 
CE phase is given by 
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where G is Newton's gravitational constant, acE is the 
CE efficiency, and Mi is the mass of the donor's compan- 
ion. Quantities with i-subscripts represent quantities at 
the onset of the CE event. 

The quantity i^bind is the difference in the donor's to- 
tal energy (gravitational potential plus internal thermal) 
between its initial and current configuration. We em- 
phasize that this represents a subtle distinction relative 
to all previous applications of this prescription. Histori- 
cally, i?bind would represent the total binding energy of 
the envelope as determined from the initial donor model. 
There is a difference between the initial envelope binding 
energy as a function of m and the donor's total binding 
energy at the corresponding total M2 along the mass 
loss sequence. This difference reflects the work done on 
overlying layers as the donor's mass shells expand under 
mass loss. This means that i^bind < for a given M2 is 
always greater than the initial binding energy at the cor- 
responding m = M2. Thus, use of the initial model over- 
estimates the amount of orbital contraction produced in 
a CE event. 

We choose a,; such that i?2,i = i?L(ai,M2,i,Mi_i), us- 
ing the expression from Eggleton (1983) to calculate Rl- 
Modulo a few complications that arise in some cases 
(e.g., systems that experience delayed dynamical insta- 
bilities) Rl{M2) is calculated during a CE event using 
a as calculated from equation Q, assuming the mass 
transfer is completely non-conservative throughout, i.e.. 
Ml = Ml,,; (see Hjellming & Taam 1991). To initiate 
the unstable mass transfer that leads to a CE event, 
(ilni?2/dlnA/2 < d\nRL/d\x\M2 initially, guaranteeing 
that R2 > Rl for some time after the onset of the CE 
event. We determine the end of the CE event to be the 
point at which R2 first becomes less than R^. 

Several points deserve to be made about this method 
for calculating CE outcomes. First, the donor's rem- 
nant mass, M2 J is an output of our model, not a pre- 
determined input. That is, our method provides a self- 
consistent determination of M2J and Of. Second, Rl is a 
function of M2, Mi^i and acE, meaning that M2J explic- 
itly depends on not only the initial binary parameters, 
but also on acE- This is in stark qualitative contrast 
with previous CE prescriptions. 

Our method is also in contrast with the existing CE 
prescriptions in terms of quantitative predictions for CE 
event outcomes. We illustrate this for the case of a 
-^2,0 = 10 M0 donor in Figure [H Here, we show M2 j 
as a function of the binary's zero-age MS mass ratio, 
q = M2.0/M1.0 and i?2,i- The systematically smaller 
M2J for smaller acE values result from a greater degree 
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Figure 2. The remnant core mass versus initial mass ratio as a 
function of binary and CE event parameters. The symbols con- 
nected by the lines denote models corresponding to a given evolu- 
tionary state (parameterized by R2,i) for which a 10 M0 donor 
star fills its Roche lobe, initiating the CE event. The solid and 
dashed-dotted line illustrates the outcomes for a CE evolution with 
aCE = 1-0 and acE = 0-5 respectively. The dotted line represents 
the core mass defined by the mass where the hydrogen mass abun- 
dance is 0.1. 

of orbital contraction. For the fiducial acE = 1-0 case, 
M2./ can vary by as much as « 20% over the q range 
producing non-merger CE events. The development of a 
deep convective envelope and steep entropy gradient near 
the core/envelope boundary on the upper RGB produces 
the transition from a sharply varying remnant core mass 
with respect to q to a more gradual decrease with in- 
creasing q seen with increasing i?2,i- 

For comparison, the dotted lines in Fig. [2] show the 
outcomes if one were to use a composition criteria to de- 
termine the core mass. We note that this value is roughly 
constant from the HG to the RGB for a 10 Mq star. A/2,/ 
determined via our method can be as much as w 30% 
greater than this core mass. 

4. COMPARISON WITH acE-PRESCRIPTION RESULTS 

To illustrate the impact of our methodology on CE 
event outcomes as compared with the standard acE- 
prescription we consider a population of binaries whose 
more massive component had an initial MS mass of 10 
M0. Some subset of this population will undergo a 
CE event as this star evolves across the HG and up its 
RGB. This initial CE event is a key step in the formation 
of many neutron star (NS) binary systems (e.g., Bhat- 
tacharya & van den Heuvel 1991; Terman et al. 1996, 
1998, Kalogera & Webbink 1998). For this comparison, 
a population of 500,000 binaries was constructed, fixing 
M2.0 = 10 M0 and considering a set of Mi_o that yield 
a uniform distribution of initial q > 1.0. The initial or- 
bital separation, oq, of this population was distributed 
evenly in Inao. For a given (g,ao) pairing, the i?2,i at 
which the 10 Mq star first filled its Roche lobe was de- 
termined. If a CE event occurred for this system we de- 
termined the resulting values of M2,/ and Of using both 
our self-consistent CE method and the standard acE- 
prescription, taking acE = 1-0 in both cases. 

The outcome of this exercise in the form of M2 j-a/ 
distributions for those systems that underwent and sur- 
vived a CE event is shown in Figure [31 The main panel 
displays the distribution of post-CE systems in this pa- 
rameter space as calculated with our method. This dis- 
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Figure 3. The post-CE population of binary systems with a 
10 Mq star tliat becomes the donor in the system's first CE event. 
The grayscale shading indicates the fraction of systems relative to 
the total surviving a CE event (in this case, 26.5% of the total 
number of binaries) as indicated by the scale on the right. In the 
two side panels, the bold lines show the corresponding histogram 
projected along each axis. The normal-weight histograms in these 
two panels show the resulting distributions when the CE outcomes 
are calculated using the standard acE-prescription. 

tributions peaks between 2.25 and 2.45 M© in M2,/ 
and below « 10 R© in af. However, both distributions 
are rather broad; the a/ distribution has a tail that ex- 
tends to significantly larger values. The comparison to 
the acE-prescription derived distributions are displayed 
in the side panels, which show the 1-D projected dis- 
tributions. Our distribution is shown by the bold-line 
histogram, the acE-prescription results by the normal- 
weight lines with the latter scaled to the total number 
of systems in our distribution. As could be easily in- 
ferred from Fig. [2l the M2J distribution in the latter 
case is nearly a J-function since the core mass of the 10 
Mq donor docs not vary significantly over its evolution 
up the RGB. The distribution in af is also significantly 
narrower using the standard prescription as compared to 
our method. 

How these differences impact the subsequent evolu- 
tionary outcomes of this binary population remains to 
be examined. However, we can provide some qualita- 
tive discussion here. Our method results in a « 20% 
higher post-CE survival rate relative to the standard pre- 
scription. This may not directly translate into a com- 
mensurate increase in the overall number of NS binaries 
produced in this channel since the core of the 10 M© 
donor will begin core He-burning and continue its subse- 
quent evolution to a core-collapse supernova. The kick 
imparted to the system during the supernova alters the 
orbital parameters and may unbind the system. Taken 
over a distribution of kick velocities and directions, bina- 
ries with larger orbital separations are more likely to be 
unbound during the supernova and have a wider distri- 
bution in post-supernova eccentricity (Kalogera 1996). 
Considering the subsequent evolution of those binaries 
that will eventually form a double NS system via both 
a second CE event and supernova (Bhattacharya & van 



den Heuvel 1991), the use of our method can be expected 
to broaden the final eccentricity distribution of this pop- 
ulation. However, our method's overall impact on the 
production rate of this (and other NS binary) popula- 
tion is unclear given the various competing processes. 
5. SUMMARY 

We have outlined a new method for calculating CE 
event outcomes that provides a self-consistent determi- 
nation of post-CE core remnant masses and orbital sep- 
arations in the framework of the acE parameterization. 
This method takes full account of the donor's structural 
response to adiabatic mass loss based on the donor's in- 
ternal entropy profile at the onset of the CE eventQ The 
essence of our method is the use of pre-calculated mass 
loss sequences to determine when, if ever, the donor is 
able to contract within its Roche radius during the course 
of the CE event. Here, the donor's final remnant mass de- 
pends on the specifics of a binary's parameters (in stark 
contrast to prior CE event prescriptions), in addition to 
assumptions about the CE efficiency. 

As an application of this method, we considered CE 
events with donors whose MS mass were 10 M©; this 
mass was chosen to highlight the significant impact our 
method has for higher mass donors that do not develop 
steep entropy gradients at their core/envelope interface 
until they are very far along their RGB. For CE events 
initiated by such objects on or before their RGB tip, ex- 
isting CE prescriptions adopt M2J to be equal to the 
donor's Hc-corc mass, which is essentially constant along 
the entire RGB. In contrast, our method predicts an M2J 
that varies by w 20% across the population of binaries 
that survive the CE event. In addition, a/ can be as 
large as « 7 times greater than the standard prescrip- 
tion. Taken over a standard distribution of initial bi- 
nary parameters, these differences result in a significantly 
broader distribution in both A/2,/ and af and an overall 
CE survival rate that is « 20% greater than the standard 
prescriptions. 

The vastly different distribution in post-CE parame- 
ters produced by this new methodology has implications 
for the multitude of binary populations of intermediate 
mass or high mass donors that experience one or more 
CE events. Our new method could be expected to pre- 
dict a broader distribution of eccentricities in double NS 
binary systems. The overall sense of impact on the pro- 
duction rate of such systems will require detailed calcu- 
lations due to several competing destruction processes. 
Future work will involve incorporating our methodology 
into such detailed considerations of this and the multiple 
other contexts in which CE outcomes play a central, if 
not uncertain, role. 
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